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PACS 81 . 16 .Rf - Nanoscale pattern formation 
PACS 68 . 37 . Nq - Low energy electron microscopy 

PACS 68 . 43 . Be - Ab initio calculations of adsorbate structure and reactions 
PACS 64 . 75 . Yz - Self-assembly 

Abstract -Spontaneously-formed equilibrium nanopatterns with long-range order are widely ob- 
served in a variety of systems, but their pronounced temperature dependence remains an imped- 
iment to maintain such patterns away from the temperature of formation. Here, we report on 
a highly-ordered stress-induced stripe pattern in a two-component, Pd-O, adsorbate monolayer 
on W(110), produced at high temperature and identically preserved at lower temperatures. The 
pattern shows a tunable period (down to 16 nm) and orientation, as predicted by a continuum 
model theory along with the surface stress and its anisotropy found in our DFT calculations. The 
control over thermal fluctuations in the stripe formation process is based on the breaking/restoring 
of ergodicity in a high-density lattice gas with long-range interactions upon turning off/on particle 
exchange with a heat bath. 
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Spontaneous formation of periodic patterns is an exam- 
ple of nature's tendency towards order. A class of such 
structures is induced by surface stress [I], and has been 
widely observed on single crystal surfaces [2HS] • It is well- 
known that the formation of these equilibrium patterns is 
driven by a competition between interactions at different 
length scales [5]. The forces in action, due to short-range 
near-neighbour and long-range dipolar interactions, are of 
the most general type resulting in very similar phenom- 
ena occuring also in magnetic [TU] and electrostatic sys- 
tems [11] . This gives a strong motivation to include such 
equilibrium phases with mesoscopic order in the discus- 
sion of making self-assembly available to future applica- 
tions [H|. 

To this date, the majority of the reports on stress- 
induced equilibrium patterns have focused on the char- 
acteristic feature shape and size. The temperature de- 
pendence of the pattern period appears as a common in- 
gredient in the various systems studied, and the general 



explanation is based on a scaling of the interaction pa- 
rameters due to thermal disorder [13] . The pronounced 
thermal fluctuations [HIE] , derived from the high mobility 
of the adatoms, are essential in reaching the thermody- 
namic equilibrium phase. However, this mobility, which 
allows the pattern to almost instantaneously equilibrate, 
is a hindrance to preserving the pattern away from its for- 
mation temperature. 

Here, we explore the possibility of controlling the ther- 
mal fluctuations by adding a second adspecies to an al- 
ready stripe-forming system. The idea is based on the 
slow dynamics in a high-density binary lattice gas, which 
leads to a tendency towards glassy behaviour. In particu- 
lar, we implement this by dosing small amounts of oxygen 
on submonolayer Pd/W(110). The oxygen-free surface is 
known to organize at high temperatures into alternating 
stripes of a dense Pd monolayer and a dilute Pd lattice- 
gas P~2] • In the following, it will be shown that the mixed 
Pd-0 layer necessitates special growth conditions in or- 
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Fig. 1: Oxygen adsorption on 0.4 ML of Pd on W(110) at constant temperature (about 1200 K). Panels (a) through (d) are 
snapshots at oxygen doses 0.54, 0.68, 1.26 and 3.51 Langmuir, respectively. The corresponding oxygen coverages are 0.08, 0.10, 
0.19 and 0.40 ML. At 5 eV electron energy, bright areas correspond to the Pd monolayer. Curved lines are W single atomic 
steps. The image size is 3.5 x 3.5 fim. 



der to form a similar stripe phase. As we will describe 
shortly, this allows to identically preserve the equilibrium 
pattern, produced at high temperature, at lower tempera- 
tures. Moreover, the addition of oxygen leads to a progres- 
sive change in the surface stress, which is reflected by the 
changes in pattern period and anisotropy. The result is a 
highly-ordered stripe pattern with a tunable period down 
to below 20 nm, and with an orientation that can be ro- 
tated by 90°. The density- functional theory (DFT) calcu- 
lations, in conjunction with continuum model results, are 
used to predict the phase diagram of stripe shape and ori- 
entation of this system as a function of Pd and oxygen cov- 
erages. This allows us to identify the equilibrium phases 
and successfully account for the main observed changes 
in stripe pattern, including the 90° rotation in the stripe 
orientation. 

Using low-energy electron microscopy (LEEM) [TlHTrj] 
we studied Pd-0 pattern formation on W(110) following 
different growth pathways and temperatures. The pro- 
gressive changes to the Pd/W(110) stripes upon oxygen 
dosing and the kinetic factors in forming the new Pd-0 
stripe pattern are revealed in the two approaches to be 
summarized, i) dosing oxygen on a partially Pd-covered 
surface at 1200 K, and ii) depositing Pd on a partially 
oxygen-covered surface at 1240 K. The latter results in 
the Pd-0 stripe phase as Pd is deposited at a tempera- 
ture at which it partially desorbs. This new stripe phase, 
with perpendicular stripe orientation relative to the Pd 
stripes on the oxygen-free surface, is identified as an equi- 
librium configuration. The conditions under which the 
Pd-0 stripes arc obtained involve effective fluctuations via 
partial desorption of Pd, which enable to reach the mini- 
mum energy configuration. 

Oxygen on W(110) forms a series of ordered phases as a 
function of coverage [T7] and has been extensively studied 
as a model system for surface diffusion and island forma- 
tion [IS]. For coverages below 0.5 ML, relevant to this 
study, it is found in (1 x2)-ordered islands. Importantly, 



oxygen has a strong tendency to laterally segregate from 
most transition metal adatoms coadsorbed on W(110) |19j . 
Furthermore, recent work has shown that oxygen induces 
a considerable stress change on tungsten [20] . These obser- 
vations point towards the possibility of finding an ordered 
two-component (oxygen and metal) adsorbate phase. 

Following this reasoning, we first dosed molecular oxy- 
gen on a partially Pd-covered W(110) surface at 1200 K, 
above the stripe disordering temperature. At this tem- 
perature, Pd shows only short-range order [13) . Upon 
adsorbing small amounts of oxygen the stripe phase ap- 
pears along [110] with a wavy morphology as shown in 
fig. [TJi. With increasing O coverage the pattern period in- 
creases (fig. QJ>) , followed by segregation of Pd into large 
2D islands surrounded by O (fig. [JJ3) . Further increase in 
O coverage (above roughly 0.3 ML) results in the shrink- 
ing of the Pd islands (fig. [TJi). This suggests that when 
the local density of O exceeds that of the (1 x2)-ordered 
layer, Pd atoms are displaced into the second layer, from 
which they readily desorb at this high temperature. It is 
important to note that at this temperature the order of 
deposition has little influence on the final result. In other 
words, depositing Pd on a partially oxygen prccovered sur- 
face below 1200 K leads to similar structures as displayed 
in fig. [EL 

In the second approach, Pd was adsorbed at 1240 K 
(about 40 degrees higher compared to that in fig. [JJ) on a 
surface precovered with less than 0.5 ML of oxygen. This 
temperature corresponds to the tail of the thermal des- 
orption peak of the first Pd adlayer [21], so that Pd ad- 
sorption is partially balanced by desorption. Under these 
conditions, a highly-ordered stripe pattern develops along 
the [001] direction as illustrated in fig.HJ Crucially, partial 
desorption was established as the necessary condition for 
forming the mesoscopic pattern. The stripe pattern con- 
sists of alternating Pd and O phases. The stripe period 
strongly depends on the amount of Pd on the surface, and 
a minimum period of about 16 nm is reached at the highest 
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Pd coverage (about 0.33 ML) allowed by the preadsorbed 

0. The new Pd-0 stripe phase differs from the pure Pd 
stripes in the orientation of the stripes ([001] vs. [110]), 

1. e. by a 90° rotation in the stripe direction, as well as 
in the improved directionality and reduced periodicity (16 
vs. 60 nm). Moreover, the pattern can be frozen upon 
cooling and thus preserved at room temperature. 




Fig. 2: Pd adsorption on 0.33 ML of O on W(110) at 1240 K. 
Panels (a) through (d) are snapshots at Pd coverages 0.18, 
0.23, 0.27 and 0.31 ML, respectively. The stripe direction cor- 
responds to the substrate [001] axis. At 13.5 eV electron en- 
ergy, dark areas correspond to regions covered with Pd. W 
single atomic steps appear as curved lines. The image size is 
1.0 yon x 1.0 fim. 



In order to elucidate whether this new phase repre- 
sents the thermodynamic equilibrium, the energetics and 
anisotropy of the Pd-0 adsorption system have been in- 
vestigated by DFT calculations, providing the parameters 
to the continuum-model arguments developed for stress- 
induced equilibrium patterns [2"2"H2"4"] . The DFT pseu- 
dopotcntial calculations were performed in a plane-wave 
basis by using the PWscf code [25] ■ The details of the 
surface stress calculations for W and O/W can be found 
in Refs. [20] and [26]. The Pd-0 boundary energy was 
calculated on a 6-layer slab consisting of a Pd-0 surface 
layer and 5 W layers. The lateral size of the slab corre- 
sponds to 10x2 substrate unit cells, of which one half is 
covered by Pd and the other half by (1x2)0. In estimat- 
ing the creation energy of a single boundary, only the first 
atoms adjacent to the Pd-0(1 x 2) boundary were allowed 
to relax (in the upper 4 layers) in order to minimize con- 
tribution from elastic interactions between neighbouring 



boundaries. 

It should be noted that the atomic structures of both Pd 
and O stripes at high temperature are a challenge for mod- 
cling. On one hand, the low-temperature ordered oxygen 
phases lose their long-range order several hundred degrees 
below the temperature at which stripes form (for example 
the (1x2) structure disorders at about 700 K [H]). How- 
ever, the short-range order persists at high temperatures, 
and influences the stripe phase profoundly. In particular, 
a local (1x2) density is the highest oxygen packing that 
can be attained in the presence of Pd islands, before Pd 
is displaced to the second and higher layers. Moreover, 
the oxygen stripe in the Pd-0 stripe phase has the (1 x 2) 
structure when cooled down to room temperature. There- 
fore, at high temperature the oxygen-covered regions can 
be described as having short-range (1x2) order. The 
effect of disorder on the oxygen surface stress is studied 
in Ref. [26], and the values obtained are reproduced in 
Tabled] 
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Table 1: Calculated surface stresses (r) and the boundary en- 
ergy (/3) between Pd and O phases on W(110). The upper 
two rows list the stress values. The lower two rows display the 
directional dependence of stress differences and boundary en- 
ergy. All stresses are given in N/m, and the boundary energy 
has units eV/A. 0(1 x 2) stress is taken from Ref. [26]. 

On the other hand, at high temperature the submono- 
layer Pd is pseudomorphic on W(110) only along [001] 
and shows a mismatch structure along [110] (i.e. one- 
dimensional pseudomorphism) |21j . Using low-energy 
electron diffraction, we observe a very similar situation 
for Pd within the Pd-0 stripes with a mismatch of ap- 
proximately 5 %. We have estimated the effect of this 
non-pseudomorphism on the surface stress by varying the 
mean lattice constant of the Pd overlayer along [110] for 
selected commensurate periodicities in DFT calculations 
using supercclls with large [110] dimensions on asymmet- 
ric slabs [57] . The calculated values for 5 % mismatch are 
displayed in the upper half of Table [T] Note that small 
variations of the mismatch (less than 1 %) due to changes 
in temperature or coverage are found to modify the cal- 
culated stresses only slightly and have no consequence in 
the following discussion. 

The results of the calculations relevant to the stripe 
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Fig. 3: (Color online) a) The phase diagram as a function of Pd and oxygen coverages generated from the calculated stresses. 
The color scale gives the Ar^jQj/Arfooi] ratio. Solid diagonal lines are contours corresponding to the critical ratio, —0.44 and 
— 1/0.44. Blue and black dashed lines trace the coverages in fig. Q] and fig. [5] respectively. Solid squares and circles mark 
the panels in each figure. Total coverage is limited to below a pseudomorphic monolayer. The Pd lattice gas is sketched for 
temperatures near Pd disordering, b) Three distinct adsorbate patterns on W(110). The patterns are qualitatively sketched 
on a diagram of stress change across the stripe boundary along the two crystallographic directions. All images have 2.0 \im 
diameter. 



formation, i.e. the directional dependence of stress dif- 
ference and boundary energy between the two phases (O 
and Pd), are summarized in the last two rows in Tabic [TJ 
For comparison, we include the stress anisotropy for the 
oxygen-free Pd stripes in the first column. An impor- 
tant finding is that the Pd-0 boundary energy is nearly 
isotropic. The values in the table clearly identify the stress 
difference as the real source of anisotropy. The clastic en- 
ergy gain is higher for a larger stress difference across the 
stripe boundary, i.e. perpendicular to the stripe direction. 
As a consequence, the expected stripe direction strikingly 
changes from [110] for oxygen-free Pd stripes to [001] for 
Pd-0 stripes on W(110). 

The calculated stress anisotropy is found to be more 
than sufficient for driving the formation of the Pd-0 
stripes along [001] according to the phase diagram (see 
fig. evaluated on the basis of a continuum-model de- 
scribed by Suo and coworkers [22TI24] and our computed 
stresses. The model predicts distinct stripe phases as the 
ratio of stress differences along the two high-symmetry 
crystallographic directions, r = At^jqj/Atjooi], is varied. 
A critical value, r c = 1 + 2v (where v is the Poisson's ra- 
tio), is found analytically by minimizing the free energy 
as a function of stripe orientation, which marks the on- 
set of regular stripes [21]. In particular for an elastically 
isotropic substrate as W(110) [2S], when r is greater than 



r c = —0.44 (or smaller than —1/0.44 for the perpendic- 
ular stripes) a linear stripe phase is expected. As can 
be seen from the results in Table [TJ this condition is sat- 
isfied both for Pd and Pd-0 stripes. In order to check 
this quantitative criterion for all coverages displayed in 
fig. the stress differences can be estimated by subtract- 
ing a weighted sum of oxygen and tungsten surface stresses 
from that of Pd. This is based on the observation that in- 
creasing Pd coverage squeezes the oxygen atoms towards a 
fully-packed (1x2) order at the highest density. The cor- 
responding stress ratios, shown as dark circles in fig. [3^, 
confirm the presence of [001] stripes for 0.33 ML O cover- 
age, regardless of the Pd coverage. 

Based on the continuum-model, a herringbone pattern 
( "0-stripes" [35]) is expected when the stress differences 
are of similar magnitude but of opposite sign along orthog- 
onal directions, with the Ar-ratio in the range —1/0.44 < 
r < —0.44. Inspecting Tabic [TJ this corresponds to an in- 
termediate situation between O-free Pd stripes and fully 
packed Pd-0 stripes. In agreement with this predic- 
tion, the herringbone pattern, already hinted by the wavy 
stripes shown in fig. [T]d, is observed in perfect form when 
the experiment is repeated, always under simultaneous Pd 
adsorption-desorption, with much lower oxygen coverage. 
The resulting pattern is displayed in fig. \$jp, together with 
the O-free Pd and Pd-0 stripe phases. These three pat- 
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Fig. 4: Stripe period as a function of Pd coverage. The fit 
function is described in the text. 



terns correspond to the three different solid-state phases 
of the overlayer in the phase diagram of fig. 

In addition to the pattern shape and orientation, the 
stripe period shows a strong coverage dependence. Using 
the expression from the continuum model [5], the period, 
D, is given by, 



D(6 Pd ,9o) = 



cos( 



(1) 



where dpd and do are the Pd and O coverages, At is 
the coverage-dependent stress difference between the fully- 
packed Pd and the partially-covered oxygen stripes, A and 
B arc the fitting parameters corresponding to the effective 
boundary width and a combination of energetic factors, 
respectively. The denominator accounts for the effect of 
stripes with unequal widths [5]. The stress difference, At, 
is found by a weighted average of the calculated values 
displayed in Table [T] as the Pd coverage is varied. The 
excellent fit to the experimental data is shown in fig. [4] 
The only fit parameters, A and B, are found to be 0.38 ± 
0.04 A and 55.6 ± 1.0 N 2 /m 2 , respectively. 

The fact that the experimental stripe period, shape and 
orientation show very good agreement with the contin- 
uum model (along with the DFT-calculated energy and 
stress parameters) confirms that the observed stripes rep- 
resent the minimum of the free energy, which is determined 
by the competition between short range (atomic bonding) 
and long range (elastic) interactions. In other words, the 
stripes are in thermodynamic equilibrium. However, as de- 
scribed in the context of the measurements, the necessary 
condition for obtaining this equilibrium phase is simulta- 
neous Pd adsorption-desorption, i.e. an open system. On 



the other hand, in the closed Pd-0 system the totally seg- 
regated state seen in fig. [TJ; is found, which is apparently 
a local minimum of the free energy. The kinetic barriers 
in this complex free-energy landscape are overcome only 
with increased fluctuations via exchange with a heat-bath 
(3D gas), allowing the system to reach the real equilibrium 
configuration, i.e. the stripe pattern. 

It is to be expected that a system with long-range inter- 
actions should face heightened kinetic barriers in its search 
for the global minimum in the configuration space. The 
basis for this statement is apparent in the formation of 
the equilibrium state. Generally, the predominant mech- 
anism in the decay towards the minimum energy phase 
is by nucleation of spatially confined droplets, which pro- 
gressively increase in size [35] . Such a possibility does not 
exist for a phase stabilized by long-range forces, which 
necessitate its formation in a coherent manner within spa- 
tially extended regions. The difference is clear in studies 
of dense and viscous lattice gas models. In particular, 
systems with short-range interactions alone do not show 
glassy behaviour, in spite of their slow-dynamics at high 
density [30]. On the other hand, long-range interactions 
in a thcrmodynamically closed system may lead to broken 
ergodicity and infinite- lifetime metastable states [3"T] . 

Importantly, it has recently been shown that the ergod- 
icity, which is broken in a closed system due to long-range 
interactions, is restored by coupling it to a heat-bath [3"2] . 
This distinction between thcrmodynamically closed (mi- 
crocanonical) and open (canonical) systems in their decay 
to the thermal equilibrium state provides a direct expla- 
nation to our observations [33] . In other words, the simul- 
taneous Pd adsorption-desorption condition is identified 
as coupling the otherwise-closed two-dimensional adsor- 
bate layer to a Pd reservoir (three-dimensional Pd gas), 
which acts as a heat-bath and allows the formation of the 
equilibrium stripe phase. 

In conclusion, we have shown the presence of a stress- 
induced Pd-0 stripe phase on W(110). The anisotropy 
and the period of the stripes depend sensitively on the 
surface stress, and can be tuned by changing the rela- 
tive coverage of the two species. Due to the presence of 
long-range interactions, the high-density lattice gas is able 
to decay into the equilibrium stripe phase only via mate- 
rial exchange with the gas phase. As a consequence of 
suppressed diffusivity, the pattern can be preserved by 
cooling to lower temperatures. This opens a new path 
towards fabricating highly-ordered periodic structures via 
self-organization at the nanometer scale. 
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